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Abstract 

We present some results concerning the large volume limit of loop 
quantum cosmology in the flat homogeneous and isotropic case. We 
derive the Wheeler-De Witt equation in this limit. Looking for the 
action from which this equation can also be obtained, we then address 
the problem of the modifications to be brought to the Friedman's 
equation and to the equation of motion of the scalar field, in the 
classical limit. 

I. Introduction . 

In recent times we have witnessed an increasing interest in the aplications 
of the ideas of loop quantum gravity (lqg) to the problems of cosmology. This 
started with a series of seminal papers by M. Bojwald [1-8] on loop quantum 
cosmology, where he obtained a number of impressive results, among them 
the possibility of removing, in a natural way, the presence of the cosmological 
singularity. 

Although importing from loop quantum gravity the techniques and ideas, 
a number of issues still remain unclear. Loop quantum cosmology (lqc) begins 
with a drastic reduction of the phase space of the full theory, resulting, among 
other things, in contrast to the situation in the full theory, in a quantization 
ambiguity in the form of a parameter which, following [9], we shall call /io, 
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playing the role of regulator. The value of fiQ canot be calculated within lqc, 
but has to be taken from the full theory, where it signals the fundamental 
discreteness of geometry, one of the most fascinating results obtained by loop 
quantum gravity. Other problems involve, for instances, the interpretation of 
the notion of quantum state and the issue of the time variable, which variable 
in lqc is replaced by the momentum variable controlling the evolution of the 
discrete quantum constraint equation (see references 10 and 11 for further 
comments on lqc). 

Another important ambiguity arises, as a result of the reduction in the 
phase space of the theory. It is well known that the classical hamiltonian 
of general relativity can be written as the sum of two terms which, in the 
homogeneous and isotropic case, become proportional to each other; the 
hamiltonian reduces itself, then, to one term. Using this property, and then 
quantizing, a discrete equation for the evolution of the system is obtained. 
This is the procedure normally used in the applications [9]. However, if we 
separately apply, to the two terms, the same quantization procedure, with 
the same connection and holonomy, we obtain instead equation (20) below 
(see also ref. [7]). Thus, it is not indifferent the order followed to derive 
the evolution equation. In addition, the two equations lead to different large 
volume limits (section III). 

An important part of the work done on lqc has turned around the study 
of the semi-classical behaviour of the solutions of the constraint equation, 
for different matter hamiltonians, in the large volume limit. This will also 
be the focus of our present work. Within the context of a homogeneous 
and isotropic universe, assumed to be dominated by a massive scalar field, 
we derive a modified form of the Wheeler-De Witt equation from the large 
volume limit of lqc, starting from equation (20). We then address the issue 
of the modifications to bring to the Friedman's equation and to the equation 
of motion of the scalar field, equations usually applied in lqc, by looking 
for the action from which, using conventional methods, we derive the above 
modified form of the Wheeler-De Witt equation. 

II. The quantum evolution equation. 
Ha. The hamiltonian for the scalar field. 

Given the importance of scalar fields in modern cosmology, we derive in 
this section the hamiltonian for the scalar field, in a flat homogeneous and 
isotropic universe, with the help of the expressions for the connection and 
holonomy given in [9]. 
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The holonomy is 



uc n . uc . 
fti = cosy + 2am— Ti, (1) 

where \x defines a length along which we integrate, not a physical length, 
and c is the real variable containing the information on the connection. The 
basis vectors, defining the eigenstates of the momentum operator p, are, in 
the bra-ket notation, 

iuc 

< c\fi >= exp (2) 

with the operator 
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P= - l —Jc (3) 



acting as 



;2 



A. I l l P l : / 1 \ 

/x e (— oo, +oo) and = 87rG. The eigenvalues of the volume operator 
are obtained from the relation V — \p\ 3 ^ 2 : 

^ = (^H) 3/a - (5) 

In the expressions for p M and V^, the constant 7 represents the Barbero- 
Immirzi parameter. 

We apply these expressions to the hamiltonian operator derived by Thie- 
mann in [12], which, in the case of a real scalar field, becomes, after the 
appropriate calculation of the traces: 



h^m = {^(villc - K^m + -fv.wmMt)- (6) 

The scalar field <fi and the momentum operator p^ = —ifid/d(f) have been 
rescaled and are dimensionless; we can check that has the correct dimen- 
sions of a mass. In the case of a scalar field of mass m, we have the potential 
W{<f>) = \m 2 (f) 2 (m in planck units). 
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We see that the operator naturally selects the ambiguity parameters, 
entering into the definition of the inverse scale factor, as I = 3/4 and j = 1/2. 
Writing, as in ref. [9], the physical state as the sum 

(*\ = Y,1>M<i*l (7) 



where represents the matter (scalar field of mass m) degrees of freedom, we 
derive the action of the hamiltonian on t/> m (</>). With a potential proportional 
to 2 , we may look for solutions in the form of linear combinations of the 
Hermite special functions. Taking the simplest case, 



we find that 
where 

and 



ipM = ^exp(-a^0 2 ) (8) 



H^M = E M<P) (9) 
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E„ = -y/A^m% (10) 



a, = ^B,/A,mfi (11) 
with the expressions for and taken from (6): 



_ 4r?v 1/2 1/2 6 

M q/9 lV+W V H-fJ.o) \ LZ ) 

and 

B, = (13) 
With these definitions we rewrite in the form 

H^^A.Pl + B^-m 2 ^. (14) 

In section III. we shall apply these expressions to integrate the Wheeler-De 
Witt equation, using numerical methods. 
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lib. The discrete evolution equation obtained from the gravitational op- 
erator. 

To obtain this equation we shall apply the holonomy operator (1) to the 
full hamiltonian operator derived by Thiemann (see Thiemann's [13] and also 
the report by A. Ashtekar and J. Lewandowski [14]), defining the regulated 
constraint operator: 

(m R£ (N) = J2(VlC e n ud (N) - 2(1 + 7 2 )Th(A0), (15) 
□ 

where the euclidean scalar constraint is 

Ch ud (N) = -ji^^e^TriCh^ - h^XKT'K^]) (16) 

and the operator corresponding to the extrinsic curvature terms 



7 AT - ~ - ~ 

Tn(N) = -jsT^ e ijk Tr ( ( /if ) ~ 1 [/if , A"] ( hf )~ l [h^\ K'}(h^y l [h^ ,V}); 

(17) 

the smearing function N will be put equal to 1 throughout and the extrinsic 
curvature operator K' expressed as a commutator [13, 14] 

K' = -^[V,C^(1)]. (18) 

Working out these expressions, the full constraint equation becomes 

(¥|(££2) + 167rGJf,)t = 0; (19) 

we shall use the expansion (7) for (\&| and, with the help of equations (15-18), 
rewrite (19) as a discrete evolution equation ( units c=l, h=l) (see [7], where 
such expression was for the first time derived): 

3 

^2^3/2^3 [(^+5W - ^+3«.)^M+4w(0) ~ 2(T^+ W - V n )^^((f)) + 

27 1 

-2(1 + 7 2 )— 9 [(V^ +9mu - V^ +7w) )ZVi_ 8w) VV+8 W) (0)- 
z * p i7 po 



-2<y M+M0 - V„)D^(<f)) + (V^ 7lM) - V^_9 W) )D M _ 8w) ^-8 W ,(0)] = 

= -IGkGH^M, (20) 

where the coefficients -D^ are complicated functions of given in the ap- 
pendix A. This equation constrains the functions Vv(^)> defined in equation 
(7), to ensure that the (\J/| belong to the space of physical states. As men- 
tioned in the introduction, the fact that this equation is different from the 
corresponding equation in reference [9] , shows that the operations of reducing 
the phase space and quantizing do not commute between them. 

We used the ordering scheme adopted by Bojwald. Had we adopted a 
symmetrized version, we could verify that Vv=o would not decouple from the 
system. 

It is known that, in the case of constraint operators, we are only interested 
in their kernel and the operators are not required to be self-adjoint. This 
seems to be an advantage in quantum gravity, an open algebra category of 
constrained system [13], as the self-adjoint requirement may lead to quantum 
anomalies. On the other hand, dropping this requirement, we are left with 
the following awkward situation. The matter hamiltonian is required from 
quantum mechanics to be self-adjoint, and indeed obeys this requirement; 
then, taking the adjoint of the evolution equation, we recover the same H^, 
but not the same H^al We are left with two distinct choices, where we 
should have one. 

However, we are here interested in the large \x behaviour of the system and 
it happens that both versions, symmetrized and non-symmetrized, lead to 
exactly the same equations; we may study this behaviour without bothering 
with this problem. To investigate the limit of large /i, it is more convenient 
to study the differential equation resulting from (20) in such a limit. 

Before proceeding, let us mention that, applying the methods developped 
in [9], section 4.3, we again find the classical limit 

< H l$>=Wj2<%y/P+--- ( 21 ) 
l pi i 

with p = 7/^/10/6 , n >> 1, where we added contributions coming from the 
two terms in equation (20). 

III. The large volume limit. 

It is important to keep present the fact that the large volume limit, large 
values of //, it is not the same as the limit, usually taken, p, — > 0. This 
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can be seen, for instances, from inspection of equation (42) in [9]. It is the 
behaviour in the large volume limit we are interested in. Working out this 
limit in terms of (/10//1), all the volume functions appearing on the l.h.s. of 
the evolution equation become proportional to ^JjT-LThjlQ times powers of 
/i . Then, adding and subtracting, after multiplication by an appropriate 
constant, the terms corresponding to the first bracket (in Vv±4mo an d Vv) °f 
equation (20), we can rewrite the l.h.s. in a form appropriate to this limit: 

1 1 _|_ ^2 g ]_ 

l.h.S. = — ) — [— (V/ X + 4 / X O^+4 W (0)- 2 \/^^( < / ) ) + V / / i - 4 / i O^-4«.(0))]- 

7 7 ft ft 



9 1 + 7 2 1 
4 7 3 ftVo 



— + 8/i Vv+8 W ,O) - 4^/1 + 4/10^+4^,(0) + 



+ 6y/ih/>n(<l>) - 4^-4/10^-4^,(0) + v / /^-8/i o ^-8w.(0))]- ( 22 ) 

What we have, within the square brackets, are the algorithms for second 
and fourth derivatives. The limit we are going to take is not /i — > 0, but the 
limit of large /i; the replacement of the discrete expressions by derivatives is 
only an approximation which becomes better and better as /i gets larger and 
larger. 

Converting the combinations of /i, /i±4/io, etc. into momentum variables, 
through equations (4) and (5), we finally find the equation replacing (20): 



f *1[/(p)v^0(p> <t>)T + ^7 2 (l + 7>65(v«P, 0))' 



w -ItoGH^M), (23) 

where the derivatives are with respect to the momentum p. The term f{p) 
is given by 



( ft7^) 2 1 
288 p 2 



m = a - ^7^-7 + • • o (24) 



and is usually put equal to 1, neglecting the small corrections. 

We see that, had we taken the limit /i — > 0, we would have recovered the 
standard form of the Wheeler-De Witt equation. This not being the case, we 
find in addition an extra term involving a fourth derivative. This term comes 
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from the two-fold application, when deriving (20), of the C eucl operator which 
appears through the definition of the extrinsic curvature operator K', in 
equation (18). Having begun with an hamiltonian formulation of a classical 
theory involving second derivatives, we end up with a quantum equation 
including a fourth derivative. 

In the vacuum case, equation (23) has analytical solutions of the form 
found by Bojowald in [7], with two of them obeying the pre-classicality con- 
dition of mild variation. When the matter term is included, we have not 
found solutions that could be put into an analytical form. We included the 
matter term in the form of a scalar field, using the procedure outlined in sec- 
tion II. a, equations (7 - 14). We solved equation (23) by numerical methods 
(see figures 1-4, with y = y/pij)) and found that we continue to have solutions, 
fig. 1 for instances, whose behaviour can be made as slow as we wish, by 
a convenient choice of initial conditions. They thus seem to fulfill the basic 
requirement for pre-classicality. 

IV. The classical equations of motion. 

We would now like to ask which classical action, using the conventional 
methods, gives us back the Wheeler-De Witt equation (23), without the term 
in (y/pip)"", and with defined by (14). (We follow J. B. Hartle's lectures, 
ref. [15]). 

The restriction to a homogeneous and isotropic geometry gives the sim- 
plest class of minisuperspace models. The line element is 

ds 2 = -N 2 (t)dt 2 + a 2 (t)dx 2 , (25) 

where now N(t) is the (arbitrary) lapse function and a(t) the dimensionless 
scale factor. We continue to take for the matter content a massive scalar 
field and define our effective classical action by the expression (representing, 
for the moment, the scalar field by tp) 

S(a, (p) = - [ dtN[-^—(-6-^-(—) 2 + 6ka) + J—(^) 2 - 

- ^B(a)m 2 pl mV, (26) 

where m is in planck units, (') = d/dt and the lagrangean is given by L = 
SS/St; from now on we put k — 0, flat case. The only non-conventional 
elements are the functions A (a) and B(a), to be defined later, and f(a), 
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given by (24), now expressed as a function of the scale a. Introducing the 
momenta 



5L 6 aa' 

na= M = ~ IQttG Nf(a) { ] 

^ = W = liv' (28) 

we find the constraint by varying L with respect to N; writing the resulting 
equation in terms of n a and ix^ we have 

^\<m - l^K - \B{a)ml^ = 0. (29) 

(Note the positions of A(a) and f(a) in (26) and (29)). 

We shall now replace (p by (f)/V6, followed by the usual replacements of 
7r a and by the operators 

** -''■Titd^ (31) 

where we took to be dimensionless. The factor ordering ambiguities, af- 
fecting the gravitational operator, will be lifted by writing it in the form 
suggested by lqc; then equation (29) becomes (|m^ = j^q) 

after dropping an overall constant factor I(j^q) 2 - 

To find the Wheeler-De Witt equation (23), all that remains is to express 
a in terms of p and introduce A(p) and B(p), the already known functions 
given by equations (12) and (13), but now expressed directly in terms of p 
(V = |p| 3 / 2 ). Remembering that a is dimensionless (a 2 is in fact the p of 
sections II and III), we introduce the relation 

~ ^ad-aad-J = ^df ] (33) 

and the definitions 
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16nGA(p) = l^fX^p) (34) 

167rG5(p) = \{-f/%^{p) (35) 
3 7 

#(p,0) = v/pV(p,0)- (36) 
With the help of these expressions we find 

f'l|^(/(p)v^(P>0) = 167rG [^^^2 - ^(p)"»W(p,0, (37) 

which is equation (23) without the the term (y/ptp)"" and with if^ given by 
(14). Thus, we may now use the action S(a, (p = 0/\/6) to derive the effective 
classical equations of motion (see Bojowald, ref. [16], for a different way of 
arriving at these modified equations). 

Again, varying L = SS/St with respect to N we get the modified Fried- 
man's equation, after replacing A(a) and B(a) by A(a) and B(a): 



1 V )2 = 8nG { l_ { 1 ( 6 )3/2](0/)2 + 1 ( 5(o^ ( 7 )3/2]m202}> 



/(a) a 3 l 2 l 3Q(l pl a) 3 A(a) 7' JVYW 2 L 2(/ pi a) 3V 6' 

; i (38) 

We see that the form of the equation is correct. With A(a) represent- 
ing, respectively, the inverse volume and the volume, apart from some con- 
stants which can be absorbed into and m, we have that a 3 A(a) and 
-B(a)/a 3 become equal tol in the large volume limit; we thus recover the 
standard form of the equation, taking also into account that f(a) as 1. 
Variation with respect to the scalar field gives 

0" - + (18(^) 3 A(a)5(a))m 2 = 0. (39) 

(See also Appendix B). Equations (38) and (39) are, with small differences, 
the equations normally used in the classical limit, when we wish to include 
corrections derived from lqc. 

That is, from the same effective classical action we can derive both, the 
modified Wheeler-De Witt equation (23) and the equations (38) and (39). 
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As for the term d 4 /dp 4 , we have not yet found a reasonable way of incorpo- 
rating into the effective action a new term corresponding to this fourth-order 
derivative. At first sight, we might think that, at least in simple situations, 
such a term would give a correction to the Friedamn's equation proportional 
to (a'/a) 4 , but a lot of care is necessary. This would come out by adding 
a term to S of the form ca(a') 4 /N 3 , c proportional to the constants multi- 
plying d 4 /dp 4 in (23) (in planck units, c is of the order of 7 2 /i(j)- The new 
momentum n a would become 

in those, and only in those, situations where the second term in the first 
bracket can be neglected, that is, when ( / -ffio) 2 l 2 l (a') 2 <C 1. Then, varying L 
with respect to N, we find that (for simplicity, we put f(a) = 1) 

167TG1 2 167TG 1 2 3 4 1 4 

7T„ -> 7T„ -c(167rG) 4 — 7T„. (41) 

12 a a 12 a a 6 4 V a 3 a V ; 

After the replacement of 7i a defined in (30), we would get, as in (23), the 
extra fourth-order derivative term y/pd 4 /dp 4 (the factor y/p later cancelling 
out) and, in the Friedman's equation, a new term proportional to a 2 (a' /a) 4 . 
Corrections of this form to the Friedman's equation do, sometimes, appear 
in the literature of brane based cosmology ([19] and references therein). 

V. Summary. 

In the present work, after deriving a modified form of the Wheeler-De 
Witt equation, as the large volume limit of the discrete evolution equation 
of loop quantum cosmology, we looked for an effective classical action from 
which this equation could be obtained, using the conventional methods. Once 
this action has been defined, we may also derive from it the associated clas- 
sical equations, the modified Friedman's equation and the scalar equation of 
motion for the massive scalar field, which we took as our matter content. 
These are the equations normally used in the applications of the loop quan- 
tum cosmology to the study of inflation. Numerical simulations were done 
with the modified Wheeler-De Witt equation, and we found that solutions 
exist that comply with the requirements of pre-classicality. We also found 
numerical solutions to the modified equation for the scalar field (Appendix 
B), having, as expected, the appropriate inflationary behaviour. 

Further applications will be left to a future work. 
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Appendix A 

The functions appearing in equation (20) are given by the following 
expressions: 

- (^+4 M0 - t / M+ 2 Mo )(V M+3w) - V^_ M )] (42) 
= --Di[(V r /1 _ 2/10 - V r M _ 4/i0 )(T/ u+/10 - V^_ 3/10 )- 

+2^2(^+6^0 - v^ +4/10 )(v r M+5/10 - V^ +/10 )— 

- - V^_ 6/i0 ) {V n - V^ 5 J], (44) 
where D 1 and _D 2 are 

-(^- V5Mo)] (45) 

D2 = [(V^+2 W) - V^)(T/ ti+ 5 Mo - V^ +W) )- 

-(^-V2 W )(W- Vmo)] (46) 

^1 — [(1^+10/xo ~~ ^m+8mo)(^+9W) ~~ ^u+5mo) — 

- (^t+8M() ~~ ^+6 W ))(^M+7mo ~~ ^m+3mo)] (47) 
-°2 — [(^M-6mo _ V r M _ 8Mo )(T/ ti _3 Mo - V^_ 7/i0 )- 
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and 



— (Vfjt-ano — V^_io MO )(V^_ 5/10 — V^_ 9w) )]. 



(48) 



Appendix B 

Let us apply a simple formal procedure, to obtain equation (39). Assume 
a regime where it is valid to use the commutators 

[4>,M = i% 

dp^/dt = ^[p h H^\ (49) 

and 

#M = -^,if*], (50) 

followed by a discretization, through the replacement of dt by a discrete step 
parameter h. We get the equations given by (14)) 

0(/i + h)= 0(/i) + hAuP+bi) (51) 

+ /i) = + hB^m 2 ^). (52) 

A little algebra shows that we can derive the following equation for 0: 

±(40* + 2h) - 20Gu + h) + fa)) = (1- ^-^ )1^ + h)- fa))- 

- A^ +h B^m 2 (p(fi) (53) 

which, again, suggests that, in this limit, we should replace the 3(d/a) term 
by the expression within the first bracket on the r.h.s. Numerical integration 
of this equation, for a large variety of initial conditions, with h ~ 4/i , shows 
that, after a period of slow increase, <f> goes through a period of very fast 
increase until it settles down into an almost constant value, a behaviour that 
is well known and is the starting point of quantum geometry inflation [17, 
18]. A simple numerical example can be seen in figure 5. 
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Fig. 1 - The graph shows the solution of equation (23), the initial condi- 
tions are y'" = y" = y' = 0, y = 1 at p = 100. We assumed a scalar field 
mass equal to 0.1 planck units. 
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Fig. 5 - Numerical integration of the discrete equation (53) for the scalar 
field. Initial conditions were 0(0) = 0, 4>'(0) = 0.1, scalar field mass=0.1 
planck units. 
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